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Abstract 



I obtain the dynamics of the continuous time quantum walk on a d-dimensional 
lattice, with periodic boundary conditions, as an appropriate limit of the dynamics 
of the discrete time quantum walk on the same lattice. This extends the main result 
of [8] which proved this limit for the case of the quantum walk on the infinite line and 
the quantum cellular automaton proposed in [1] . By highlighting the main features 
\^ ■ of the limiting procedure, I then extend it to general graphs. For a given discrete 

0^ . time quantum walk on a general graph, I single out the type of continuous dynamics 

(Hamiltonians) that can be obtained as a limit of the discrete time dynamics. 

1 Generalities on quantum walks; continuous and dis- 



crete time 



Consider a graph G := { V, E} with a set of vertices V of cardinality N and a set of edges 
E. We assume G to be undirected and without self-loops. In the interval of time At a 
certain fraction 7 At of a quantity pj leaves the location of the j-th vertex to move to a 
neighboring vertex k. The quantity pj, j = 1, . . . , N, at time t + At is given by 1 

Pj {t + At) = Pj (t) - deg (jhAt Pj (t) + l A tPk(t)- (1) 

(k,j)eE 

Define p the iV— vector whose components are the p/s quantities, and L the Laplacian 
matrix defined by Ljj = — deg (j), = 1 if (j,k) G E and = otherwise. L is 
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1 The rate 7 may depend on time and we could have allowed it to depend on the location j. In the 
latter case though we cannot choose a Laplacian matrix L (see later) to be symmetric. 
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symmetric and, except for the elements on the diagonal, coincides with the adjacency 
matrix A, that is, Ajj, := Ljj,, when j ^ k and Ajj, = if j = k. Taking the limit At — > 
in (1) one obtains the differential equation for the classical continuous-time random walk 
(CTRW). 

j t P = lLp. (2) 

The reason to call this model a random walk is that the p 3 -(t)'s may represent the proba- 
bility of a walker to be in position j at time t. A continuous-time quantum random walk 
(CTQW) (see [4]) is obtained by quantizing equation (2). One replaces p with a complex 
iV-vector representing the state of a quantum system, while L is taken as the Hamiltonian 
determining the evolution according to the Schrodinger equation 

iij; = 1 Lifj, (3) 

with \ifjj\ 2 representing the probability of the quantum system being in the basis state 
\j). This definition can be generalized by using an Hamiltonian different from L which 
respects the topology of the graph G, that is elements different from zero correspond to 
edges in the graph. For example one can use the adjacency matrix A instead of L in (3). 2 
The classical discrete-time random walk (DTRW) is obtained by discretizing equation 
(2) and therefore it is given by equation (1). The CTRW is obtained from the DTRW 
by taking the limit At — > in equation (1) as seen above. To define the discrete time 
quantum walk (DTQW), we would like to write an equation of the form 

ip(t+i) = um, (4) 

with if) the state of a quantum system and U a unitary operator representing a (closed) 
quantum evolution. We also would like to have a U which respects the structure of the 
underlying graph G, i.e., Ujk, with j ^ k is different from zero if and only if there exists 
an edge (j, k) in E. Unfortunately, only special graphs have the property that such a 
unitary matrix exists [7] . In order to give a definition which is suitable for any graph one 
proceeds as follows. 

Let V be the Hilbert space spanned by the orthonormal states with j e V. 

Let S be the subspace of V <8> V spanned by \j, k) with (j, k) G E. The basis state \j, k) 
represents the state of a walker which is currently in vertex j and is moving to vertex k, 
i.e., j and k are the present and future location of the walker, respectively. On the space 
S, the evolution U of the DTQW is of the form U = WC. The operation C, called coin 
tossing, is of the form 

6=EL7>0'I®Q;, (5) 
jev 

where Qj is a unitary transformation on V, depending on j. The subspace spanned 
by the states corresponding to the neighboring vertices of j (and therefore its orthogonal 

2 In this case, if G is a regular graph, i.e., deg(j) is the same for every j in V, the corresponding 
dynamics would differ from the ones of (3) only by a physically unimportant phase factor. 
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complement) is invariant under Qj. The operation W is any unitary which transforms the 
elements \j, k) as W\j, k) = \k,r), for some \k,r) G £, i.e, it moves the future state in the 
present state position. One possibility is the swap operation, defined by W\j, k) = \k,j), 
since \j, k) G £ <-> \k, j) G £. 

If the graph G is regular, one can give a definition which makes the coin's role more 
transparent. We call this quantum walk the coined DTQW. It is defined as follows. 
Let m be the degree of G and consider a (coin) space C spanned by orthogonal states 
{|ci), . . . , \c m )} each representing the result of a coin tossing. Denote by nj(c k ), j = 
1, . . . , N, k = 1, . . . , m, an element in V if a coin result c k induces a transition from j 
to it. The coined DTQW evolves on C <E> V as R = SC <g> 1, where C is a unitary (coin 
tossing) operation on C, 1 is the identity on the walker space V and S is a controlled shift 
defined by 

S\c k ,j) = \c k ,n j (c k )). (6) 

The coined DTQW, defined in the case of a regular graph, is a special case of the more 
general DTQW defined above, when the coin operations Qj in (5) are essentially indepen- 
dent of j and W takes a special form. To be more specific, let C in the coined DTQW be 
defined by 

m 

c \ c k) ■= J2 a ik\ c i)- (7) 
1=1 

Then choose Qj in (5) as 

m 

Qj\nj(c k )) = J2 a ik\ n j( c i)) ( 8 ) 
i=i 

(notice the coefficients cti k are independent of j). The transformation W of the DTQW 
is chosen as 

W\j,nj(c k )) = \nj(c k ),n nj{ck) (c k )). (9) 

With these choices, the dynamics of the DTQW on 8 are equivalent to the dynamics of 
the coined DTQW on C <E> V. In fact, using the isomorphism x '■ ^ ®V ^> £ defined by 

X\ck,j) = \j,rij(c k )), (10) 
we have that x>S(C <8> 1) = WCx- In fact, a straightforward computation shows 

m 

xS(C®l)\c k ,j) = ^2ai k \n j (ci),n nj(C[) (ci)} = WCx\c k ,j). (11) 
i=i 

Having defined the DTQW and the CTQW the question remains on whether the 
CTQW can be obtained as a limit of the DTQW on the same graph. This is a fundamental 
question which was posed in [6] and has received much attention recently (see, e.g., [2], [8] 
and references therein) . It is motivated, among other things, by the fact that CTQW and 
DTQW on the same graphs have showed similar behavior in several applications. Since 
the DTQW evolves on an higher dimensional space as compared with the CTQW, the 
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procedure to obtain the CTQW from the DTQW should involve not only an appropriate 
limit but also some sort of a projection of the dynamics onto a lower dimensional space. 
This work was done in [8] for the quantum walk on the infinite line. In this paper, I 
extend and somewhat simplify the calculation in [8] to prove the same limit for the coined 
DTQW on a d-dimensional lattice with periodic boundary conditions. In particular for 
d — 1, I obtain the CTQW from the DTQW for the cycle. I then observe that this 
generalization is only an example of a general procedure that I then extend to general 
graphs. I formalize this procedure by describing the set of Hamiltonian dynamics for the 
CTQW which can be obtained by a limit of the dynamics of the DTQW (cf. Theorem 
1). In the following, I shall consider only coined DTQW. 

The paper is organized as follows. In Section 2 I describe in detail the quantum walk 
on a rf-lattice. The d-lattice is the Cartesian product of d cycle graphs. If the transition 
from discrete to continuous is obtained for each factor of a product graph then it carries 
over to the whole graph (cf. end of Section 2 and Remark 2.1). Therefore, it is enough 
to restrict ourselves to the cycle. I prove the transition from DTQW to CTQW for the 
cycle in Section 3. The calculation in this section can be carried over to other quantum 
walks. Moreover, even for the case of the cycle alternative limiting procedures to obtain 
the CTQW from the DTQW can be devised as discussed in Remark 4.1. In Section 4, I 
extend the procedure to general coined quantum walks. I characterize in Theorem 1 the 
set of Hamiltonians on the space £ <g> V whose dynamics can be obtained as a limit of the 
one of the DTQW. I give an example in section 5 and conclude in section 6. 

2 Quantum walks on d- dimensional lattices 

Consider a coined DTQW on a d-dimensional lattice whose set of vertices V is given 
by V = {0,1,..., N - l} d . If V := span {|0>, |1), . . . , \N - 1)}, is the Hilbert space 
associated with {0, 1, . . . , N — 1}, the walker Hilbert space for this graph is V = ® d V, 
which is spanned by the orthonormal vectors \ji, ■ ■ ■ , jd) , ji, ■ ■ ■ ,jd = 0, 1, , . . . , N — 1, 
which represent vertices of the graph with coordinates ji, . . . ,jd- In the graph G, the 
vertex labeled by . . . ,jd) has 2d neighbors each differing by . . . ,jd) by only one 
coordinate and with (Hamming) distance 1, i.e., of the form (ji, . . . , 1, . . . ,jd), where 
the ±1 operation has to be intended mod N. Let G be the graph representing the cycle 
with N nodes. G is the product of d copies of 67. Its adjacency matrix is 

A:=Y,A®, (12) 
i=i 

where A® is the tensor product of d, N x N identity matrices except in the Z-th position 
which is occupied by the adjacency matrix of the cycle A. 3 A is given by A = F + F T , 

3 There are several definitions of product of d graphs. Here we intend the Cartesian product of 
d graphs Gi = (Vi,Ei), I — l,...,d which is defined as the graph G = (V,E) having the set of 
vertexes V equal to the Cartesian product V := V\ x V2 x . . . x Vd and set of edges E such that 
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where the circulant matrix F is defined as 
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(13) 
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i o/ 



The Laplacian is given by 



L :— — deg (G)l + A = -2dl + A, 



(14) 



and the solution of the corresponding Schrodinger equation is tp{t) = X(t)if)(0), with 
X{t) — Xi(t) <S> • ■ • <S> Xd(t), where Xi is the solution of the Schrodinger operator equation 
corresponding to the l-th graph (cycle) 



iX^iLXu X,(0) = 1, 
where L = —21 + A is the Laplacian corresponding to a single cycle. 



(15) 



To construct a coined DTQW on the <i-dimensional lattice, we have to introduce a coin 
space C which is spanned by 2d basis vectors, each corresponding to a different coin result. 
We choose the basis so that the vectors corresponding to the same degree of freedom are 
placed one after the other. For example, for a 2-dimensional lattice, C is spanned by the 
ordered basis {| — >), | <—), | |), | I)}, which induce a right, left, up and down motion, 
respectively. The dynamics of the coined DTQW on C <g> V is given by SC <E> l N d xN d 
where C is a coin transformation on C and S is the controlled shift. S is a 2dN d x 2dN d 
block diagonal matrix with d, 2N d x 2iV (i -dimensional, blocks each corresponding to one 
degree of freedom in the coin space. The Z-th block is given by 



(F® 
I F« T 



(16) 



where (F^ T ) is the tensor product of d, N x N identity matrices, except in the Z-th 
position which is occupied by F (F T ) in (13). It represents a forward (backward) motion 



((u\,U2, ■ ■ ■ , Ud), (vi, V2, ■ ■ ■ , Vd)) € E if and only if all Uj's are equal to the corresponding v/s except 
for exactly one pair wj ^ which is such that (uj,Vj) e Ej. The adjacency matrix of such a graph 

is A = where A^> where A^> is the tensor product of d matrices all equal to the identity, 

with the one in a generic position r of dimension \V r \, except for the matrix in position I which is equal 
to the adjacency matrix of the l-th graph. A quick way to see that this is the case is to recall that a 
definition of adjacency matrix A given by ejAek — 1 if and only if (j, k) is an edge and = otherwise. 
Here and are the standard basis vectors corresponding the vertexes j and k. A quick calculation of 
(ji, . . . , jd\A\ki, . . . , kd) shows that this is the case for the above defined matrix. If the graphs G\,...,Gd 
are all regular, so is the product graph and its degree is the sum of the degrees of the graphs G\, . . . , Gd- 
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in the degree of freedom identified by I (for example I = 1 (right or left motion), I = 2 
(upward and downward motion). 

The matrix C is, in principle, a general matrix in U(2d) (i.e., 2d x 2d unitary), which 
depends therefore on (2d) 2 real parameters. It is the only element in the dynamics which 
couples the various degree of freedom of the motion of the DTQW. We would like to 
obtain the dynamics of the CTQW as an appropriate limit involving the parameters in 
this matrix and going from discrete to continuous time. Since the dynamics of the CTQW 
given by (15) is completely decoupled, it is reasonable to restrict our attention to coin 
transformations which do not couple the various degrees of freedom of the coin. The 
corresponding matrix C has a block diagonal structure with d, 2 x 2 blocks belonging to 
U(2). With this restriction, the dynamics of the DTQW also are completely decoupled 
and we can restrict our attention to only one factor in the product graph, i.e., a cycle. 
Therefore, in the following, we shall restrict ourselves to a DTQW on a cycle with N 
vertexes, whose dynamics is given by S(C <8> Inxn) with C in SU{2) and 

Mo £)■ ™ 

Remark 2.1 In what we have said above, there is no reason to consider the product of 
cycles with equal number of nodes N, other than notational convenience. More generally, 
for product graphs, the adjacency matrix always has the form of a sum of adjacency ma- 
trices as in (12) and the dynamics of the continuous quantum walks are always decoupled. 
Therefore if we have a limiting procedure to obtain the dynamics of the CTQW from the 
DTQW for each factor graph, we can can carry this over to the product graph by restrict- 
ing the dynamics of the coin space in a way that the dynamics on the different factor 
graphs are decoupled. The same holds for the quantum walk where the Hamiltonian is 
given by the Laplacian since this only differs from the one corresponding to the adjacency 
matrix by a physically unimportant multiple of the identity. The important fact is that 
the Hamiltonian for the product graph has the form (12) (i.e., sum of matrices which 
are products of the identity on the various spaces except for a local Hamiltonian). There 
are however Hamiltonians that respect the structure of the product graph but cannot be 
written this way. For these Hamiltonians the graph has to be considered as a whole. 

3 Continuous time quantum walk as limit of discrete 
time quantum walks on the cycle 

Consider the coined DTQW on the cycle and write the coin operation C G SU (2) as 

C = Re iDx , (18) 

with R := ~^ and D := ( ^ With x small, U 2 (x) = SC <g> ISC <g> 1, 

represents a small perturbation to a transformation U 2 (0) = — 1 2 <8> ljv (cf. (17) using 
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the fact that FF T = 1 N ) which does not modify the state of the system (except for an 
unimportant phase factor). Write e lDx = 1 2 + iDx + 0(x 2 ). We write U 2 (x) as 4 

U 2 {x) = S (R(l 2 + iDx + 0{x 2 )) <g> 1 N ) S (r(1 2 + iDx + 0{x 2 )) <g> 1 N ) (19) 

= (S(R <g> l)) 2 + i ((S(R <g> l)) 2 D <g> 1 + S(R <g> 1)D <g> 15(i2 <g> 1)) a; + 0(:r 2 ). 

We have already said that (S(R<S> l)) 2 = — 12 <E> ljv- On the other hand, a direct calcu- 
lation shows that 

(S(R® 1)) 2 D® 1 + S(R® 1)D® 1S(R® 1) = ( 1+ ° F 2T 1+ F2 )- ( 20 ) 
Therefore, we have 

l^) = -(l-*( 1+ V 1+ F2 )x + O(x 2 )). (21) 
If H is the Hamiltonian 

1+ o F ^ (22) 

then 

£/ 2 ( x ) = -e-^ + O^ 2 ), (23) 

i.e., for small perturbations x two iterations of the DTQW give an evolution which (ex- 
cept for the unimportant phase factor —1) corresponds to the continuous evolution by 
a Hamiltonian H in (22) in the same time x. One obtains an evolution e~ nHt over an 
arbitrary time t by applying U 2 (x) an infinite number of times and, at the same time, 
letting x — > 0. More specifically, setting \ = and neglecting the, physically irrelevant 
factor —1 in (23), we have 5 

lim U T (x) = lim (e~ iHx + 0(x 2 )) * = e~ ilHt . (25) 

x^O x^O V ' 

Therefore, over the interval [0, t] the system evolves according to the Schrodinger equation 
iip = ryHtp where H is defined in (22). Defining ip := [^,ipJ / ] T : with ip R and ip L both 
A/-dimensional, we have 

j>R = -il{l + F a )ip L , (26) 



Sometimes, when it is clear from the context, we omit the dimension of the identity matrix 1. 
5 This limit can be obtained from general properties of the logarithms of matrices (cf. [5] section 6.5). 
Set := m = ^. The limit in (25) becomes 

lim ( e - iH %+0{±)) m = lim ( 1 - iH* + 0{ \ )) ™ = lim e^-^°^T = 

ra— >oo y ui ) m— >oo \ m 771 J m— >oo 

lim e - m {U iH ^- ^))+^( iH ^- ^)T+-) = e -^ Ht . (24) 

m^oo 

Here we have used the series expansion of the principal logarithm of a matrix log(l — A) = — J2kLi jA k - 
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i> L = -i 7 (l + F 2T )ifj R . (27) 

There are special (linear) combinations of ip R and ipi which evolve according to i$ = 
±^Afy where A = F + F T is the adjacency matrix of the cycle graph. In particular, 
define 

*i+ := Vfc + FipL- (28) 

Using (26) and (27) 

*i+ = ~il ((F T F + F 2 )^ L + F(l + F 2T )^ R ) = (29) 

-ry (F T (ij R + Fij L ) + F{^ R + Fifj L )) = -i 1 A^ 1+ . 
With analogous calculations, after defining 

:= + F t ^j r , (30) 

:= ip R - FtpL, (31) 
tt 2 - := ipL - F T ip R , (32) 

one finds 

i*i, 2 ± = ±7^i, 2 ±. (33) 
From these equations, it follows that, defining 

$i,2± := *i )2 ±, (34) 

we have 

z$i, 2 ± = ±7L$ 1)2± , (35) 

where L is the Laplacian L = —21 + A. These are the dynamics (forward or backward in 
time) of the corresponding CTQW. 

The dynamics of [iPr,i/)l] t are obtained as the average of dynamics of the CTQW 
corresponding to the adjacency matrix, i.e., (from (28), (30), (31), (32)) 



(36) 



4 Extension to general graphs 

The calculation presented in the previous section follows the main steps of [8] which we 
have adapted to the case of the cycle. We have however modified this calculation in 
several respects. In particular, we have omitted the Fourier transformation (and anti- 
transformation) step and replaced the calculation (8) of [8] which is based on the algebra 
of Pauli matrices with a Taylor expansion in a perturbation parameter x. In this form, the 
treatment highlights the main ideas of the process to obtain the dynamics the CTQW as 
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an appropriate limit of the coined DTQW and therefore can be generalized. In fact, we will 
show in this section how this limit can be extended to quantum walks on general graphs 
and-or CTQW with Hamiltonian different from the Laplacian or the adjacency matrix. 
The main idea of the calculation in the previous section is to take a reference trajectory of 
the coined DTQW which takes the walk back to its original position (up to an overall phase 
factor). 6 Then one perturbs such a trajectory by slightly modifying the coin operation 
at each step. Let x be a parameter which measures the magnitude of this perturbation. 
The resulting trajectory will agree with e~ tHx for a certain simulable Hamiltonian H up 
to higher order terms 0(x 1+6 ), (5 > 0). Then one repeats this trajectory a number of 
times which increases as x — > while letting the perturbation go to zero as in (24). The 
result is an evolution of the type e~' ljHt for some scalar 7. The Hamiltonian H acts on 
a c x N space where c is the dimension of the coin space and TV is the dimension of the 
walker space (i.e., the number of vertices of the graph) while we would like to obtain an 
evolution according to an Hamiltonian H acting on N- dimensional space. However, if the 
simulable Hamiltonian H has N eigenvalues coinciding with the ones of H, by a change 
of coordinates we can isolate a subspace of C <g> H where the evolution coincides with the 
one determined by H. This is meaning of the calculations in (28)-(36). 7 

Remark 4.1 We remark that the above procedure can be used not only to obtain the 
CTQW as a limit of the coined DTQW for new graphs (as we shall see in the following) 
but also to obtain alternative procedures for the case of the cycle. For example, consider 
the reference trajectory S N = 1 with S given in (17) and perturb it as 

gN ^ Se Ex 1S N-l e Ex 1; (37) 

with E :— ^ ^ q^J. A straightforward calculation shows that we have 

Se Ex <g) lS N ~ l e Ex ® 1 = 1 - iHx + 0(x 2 ) = e~ iHx + 0(x 2 ), (38) 

where H is the same as in (22). Therefore the treatment then goes as after formula (23). 

The question at this point is the characterization of the set of simulable Hamiltonians 
for general coined DTQW. To this purpose we give a more precise definition of simulable 
Hamiltonian and then describe this set in Theorem 1 below. 

Definition 4.2 An Hamiltonian H on the Hilbert space C <8> V (with c = dimC and N = 
dim V) is called simulable if there exists a parameter x and a sequence of transformations 
SCje E J x <g> 1, with Cj G U(c), Ej e u(c), j = 1, . . . ,m, such that 

m 

J] SC je E ^ x) 01 = e~ iHx + 0(x 1+s ), (39) 

3=1 

6 Note that such a trajectory always exists. As the matrix S in (6) is a permutation matrix S r = 1 
with r the order of S. That is: by taking r steps with the coin operation C equal to the identity we come 
back to the original position of the walk. 

7 In this case, in fact, there are 4 possible subspaces and, at the limit, the dynamics of the DTQW 
gives 4 copies of the dynamics of CTQW, 2 forward in time and 2 backward in time. 
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for some 5 > 0, strictly increasing smooth functions, fj, fj(0) = and x G [0, e), for some 
e > . 

If H is simulable, then one can use a sequence of steps of the DTQW and perform 
a limit procedure as in (24) to obtain the dynamics corresponding to H. In particular, 
using (39), one obtains 




Theorem 1 Let r be the order of the matrix S of a coined DTQW, and consider the set 

T = {u(c) <g> 1 N , Su(c) <g> l N S r -\ S 2 u(c) <g> l N S r ' 2 , S^uic) <g> 1 N S}. (41) 

Then, the set of simulable Hamiltonian is (modulo i) the Lie algebra generated by T . 

We denote by 'Sim' the set of simulable Hamiltonians and by C the Lie algebra gen- 
erated by JF. 8 The statement of the Theorem says that 

i Sim = C (42) 

Proof. First, we show that i Sim C C Assume H G Sim so that (39) holds. For every 
x in [0, e) the left hand side belongs to the Lie group corresponding to the Lie algebra C, 
which we denote by e c . To see this, we first observe that e c has the property 
Property (a): 

X G e c <=> S k xs- k G e c , 

for every integer k. 

This property holds for the Lie algebra £ (i.e., with e c replaced by C) and carries over to 
the Lie group e c . 9 

Notice that the left hand side of (39) (as well as the right hand side) is equal to the 
identity when x — 0. If m = 1 in (39) then the left hand side is a curve in e c for every 
x G [0,e), T = T(x), with T(0) = 1. This is true even if m > 2. To see this, we 
proceed by induction on m, using Property (a). From YVJLi ® ljv = 1 c at, we obtain 
SCi (8)1 = 11?=™ (cj<g)lSt), and therefore we have 

n sc je E ^ x) ® i = n (cj ® i^ f ) e si/iw ® i n ^c je ^^ w ® 1 : = L m (x). (43) 

8 Note this coincides with the dynamical Lie algebra studied in [3] which determines the set of reachable 
states for a DTQW on the cycle seen as a control system. 

9 To prove that the property holds for C just observe that it holds for the set T in (41) and it is 
preserved under linear combination and Lie bracket of two elements of the Lie algebra. 
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Define L±(x) = e El ^ x ^ <g> 1 which is in e c . We have 

L m+1 (x) = Cl +1 <g> lS^L m (x)SC m+1 <g> i e ^+i/ m+1 (-) ® 1. (44) 

S^L m (x)S G e c because of Property (a) and taking the product with element in e c we 
obtain an element in e c because of the group property. In conclusion the left hand side 
of (39) is a smooth curve T(x) in e £ , with T(0) = 1. From 

T(x) = e~ iHx + 0(x 1+5 ), (45) 

taking the derivative with respect to x at x = we obtain —iH = T(0), i.e., —iH is an 
element of the Lie algebra C. This conclude the proof that z Sim C C. 
To prove that C C i Sim we prove two facts 

1. 

T C z Sim 

2. z Sim is a (real) Lie algebra. 10 

From these two facts, since, by definition, C is the smallest subalgebra of u(cN) containing 
J 7 , it follows that C C i Sim . The proof of 1. is almost obvious. If —iH G T then 
e -iHx _ gk e Lx igr-k f Qr gome — 1 ^ r ^ anc [ this can be written in the form (39) by 

choosing m — r, Cj — 1 for every j, i?j = L for j — k and £j = otherwise, and fj(x) = x, 
for every j. To prove 2., we need to prove that: (a) iH G z Sim <^ — zi7 G z Sim ; (b) 
iH G z Sim =^ azif G z Sim , for each positive a; (c) iH 1: iH 2 G z Sim ^> z'(i^i + #2) G 
z Sim ; (d) iH 1} iH 2 G z Sim ^> [zi^x,?^] G z Sim . That is, (a), (b) and (c) show that 
z Sim is a vector space and (d) that is closed under commutation relation. Rewrite (39) 
as 

e' lHx = T(x) + 0(x 1+5 ), (46) 

(cf. (45)), i.e., by replacing the product in (39) by the symbol T(x), and notice that T _1 (x) 
can be also written as a product of the form in (39). By straightforward manipulations 
of (46) we obtain T~ 1 (x) = e iHx + T-\x)0(x l+s )e iHx , that is e iHx = T~\x) + 0(x iu ), 
which shows (a). To show (b), notice that if T(x) is an admissible product so is T(ax), 
for a > 0, so by replacing x with ax in (46), we obtain the result. To prove (c), denote 
by Ti(x) and T 2 (x) the product in (46) corresponding to H\ and if 2, respectively and 
Si and 62, the corresponding 5's in (46). From the exponential formula e~ lHlX e~ lH2X = 
e -i(H 1+ H 2 )x + o( x 2 ), we obtain 

e - l{Hl+ H 2 )x + Q ( x 2j = (j,^ + (x 1+Sl )) (T 2 (x) + 0(x 1+S2 )) , 

which gives 

e -i(H 1+ H 2 )z = Tl (x)T 2 (x) + 0(x 1+min(1 ' <5l ' <52) ). 

10 All the Lie algebras we are considering are subalgebras of the real Lie algebra u(n) (for appropriate 
n). 
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which shows our claim since T{T 2 is also an admissible product. Finally, to prove (d) 
we use the exponential formula ^r lH ^- lH ^ t _|_(9^ 3 ) — e -'iH 1 t e -tH 2 t e %H 1 t e iH 2 t ^ se6) e g ^ [5^ 

Section 6.5.). Using this and the notation described above, call Oi 2 the O function for 
Hi and H 2 , respectively. We have 

e[ -iH lt -iH a ]1* + = 

(Tx(t) + 1 ) (T 2 (t) + 2 ) (Tr\t) - Tf^O^') (T 2 ~ l (t) - T^{t)0 2 e iH ^) . 

Expanding the right hand side and omitting explicit terms that are clearly 0(t a ) for a > 2 
since they contain products of two O functions, we obtain 

e [-iff! -iff 2 ]t 2 + o(f 3 ) = 0(t a ) +TiT 2 Tf 1 T 2 " 1 -T^T^T^O^ 1 " 21 -T 1 T 2 T^ 1 1 e iHlt T 2 1 + 

(47) 

T x 2 T^T 2 x + OxTaTf 1 ^" 1 . 

Developing in a Taylor series the functions that multiply the O functions on the right 
hand side of this expression gives cancelations which show that only terms of the type 
0(t 2+Sl ) and 0(t 2+&2 ) possibly remain. In conclusion, we have 

e l~iH 1 ,-iH 2 ]t 2 = Tl (t)T 2 (t)T{\t)T 2 \t) + 0{f), (48) 

with (3 > 2 and by setting x = t 2 , we obtain 

e [-iH u -iH 2 ]x = Tl {^)T 2 {^)T{\^)T 2 -\^) + O(xl), (49) 

which show the result because if T(x) is an admissible product so is T(^fx). This concludes 
the proof of the theorem. □ 



5 Example 



I consider now an example of a coined DTQW whose associated CTQW can be obtained 
as an appropriate limit. This example is not a DTQW on cycle or ci-lattice and shows 
the generality of the method described above. In particular, consider the graph in Figure 
1 which has a three dimensional coin space C (spanned by the coin tossing results 1, 2, 3) 
and a 4-dimensional walker space spanned by the positions A, B, C, D, in that order. 
The associated matrix S is given by S = diag (Si, S 2 , S3), with 



Si 



[0 1 0\ 

1 

10 

\o 1 oj 



fO 1 0\ 

10 

1 

\o 1 0/ 



s. 



f° 
















1 








1 








U 








oy 



(50) 



The Lie algebra of simulable Hamiltonians i Sim is generated by 

T = {u(3) ® 1 4 , Su(3) ® US}, 



(51) 
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Figure 1: Graph for the discrete time quantum walk in the example. The labels on the 
edges indicate the corresponding coin toss result. In this example the same coin toss 
result induces motion in both directions. 



since S has order 2. It contains, in particular, the element diag (— 3i, i, 2i) <g> 1 4 . The 
adjacency matrix A := Si + S 2 + S 3 has one eigenvalue equal to 3 with multiplicity 1 
and one eigenvalue —1 with multiplicity 3. Therefore, we can single out a 4-dimensional 
vector space of C £g> V of the DTQW where the dynamics coincide with the one of the 
CTQW corresponding to the adjacency matrix. We can, in fact, do better by studying 
more closely the structure of the Lie algebra generated by T in (51). S\, S 2 and S3 can be 
simultaneously diagonalized, and moreover we have SiS 2 = S 3 , S 2 S 3 = Si and S3 Si = S 2 
and Sf — S 2 — Sf = 1 4 . The set T can be written as 

T = {u{3) (g 1 4 , E 12 <g> S 3 , E 13 <g> S 2 , E 23 ® Si, } (52) 

where Ej k represents the set of all the matrices in su(3) with all the entries equal to zero 
except for the j, k and j, k entries. From this, one finds that 

C= span {u(3) ®U,su(3) ® Si,su(3) ® S 2 ,su(3) ® S 3 .} (53) 

It is in fact clear that this is a Lie algebra containing T . Moreover by calculating the set 
0£l o ad k su( ^Ej k , for a given jk it is clear that this is all of su(3), since it is a non-empty 
ideal in su(3), which is a simple Lie algebra. Therefore, we obtain with repeated Lie brack- 
ets all the matrices in su(3) (g Si^^. At this point, it is convenient to make a change of 
coordinates to simultaneously diagonalize Si, S 2 and S 3 and to change the order of the ten- 
sor product in (53). Therefore the matrices in C are linear combinations of matrices of the 
form Li := diag (A, A, A, A), L 2 := diag {B,B,-B,-B), L 3 := diag {-C,C,-C,C), 
L 4 := diag (— D, D, D, —D), with A E u(3), and B,C,D E su(3). We can choose a linear 
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combination of £1,2,3,4 to make the resulting matrix have eigenvalues ±3i with multiplicity 
1 and ±i with multiplicity 3 and with multiplicity 4, so as to obtain 2 copies of the 
CTQW (one backward and one forward in time) and no dynamics on appropriate four 
dimensional subspaces of C <8> V. 

6 Conclusions 

I have obtained the dynamics of the continuous time quantum walk on the cycle as an 
appropriate limit of the coined discrete time quantum walk on the same graph. This result 
can be extended directly to the case of (i-lattices using the fact that they are products of 
cycles. The main ideas of the procedure can be extended to general graphs. One can use 
the discrete time quantum walks to obtain the dynamics corresponding to a given Hamil- 
tonian on the whole (coin+walker) space and then obtain the dynamics of the continuous 
time quantum walk by restricting oneself to a subspace of the whole (coin+walker) space. 
I have characterized the set of Hamiltonians for which this is possible, a set that I called 
of simulable Hamiltonians. The results of this paper reduce the problem of obtaining the 
continuous time quantum walk from the discrete time to the problem to a study of the 
set of simulable Hamiltonians. This set has the structure of a Lie algebra. As this set is 
quite rich it is reasonable to expect that every (or at least most) continuous time quantum 
walks can be obtained by a limiting procedure and restricting to an appropriate subspace 
the dynamics of the associated discrete time quantum walk. 
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